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FOREWORD 


This  report  was  prepared  under  Air  Force  Contract  Number 
AF  08(635 )-2155,  "Study  of  Target  Penetration  Prediction  By  High 
Speed  and  Ultra  High  Speed  Ballistic  Impact".  Work  was  administered 
under  the  direction  of  APGC  (PGTWR),  Eglin  Air  Force  Base,  Florida, 
with  Mr.  A.  G.  Bilek  as  Project  Engineer. 


Catalog  cards  may  be  found  in  the  back  of  this  document. 


ABSTRACT 


A  statistical  analysis  of  the  penetration  depth  in  semi-infinite 
targets  divided  by  the  diameter  of  the  projectile  (Pc/Dp)  on  1272  ex¬ 
perimental  shots  has  been  completed.  This  data  was  split  into  985  low 
velocity  shots  and  297  high  velocity  shots  according  to  the  bulk  wave 
velocity  in  the  target  material.  Separate  analyses  of  the  two  groups 
show  interesting  relationships  with  existing  theoretical  and  empirical 
equations. 

A  semi-rational  penetration  expression  has  been  developed  from  a 
work-energy  consideration  which  suggests  that  the  nonrecoverable  target 
compression  and  shear  strain  energies  may  account  for  most  of  the  kinetic 
energy  of  the  projectile.  Judging  from  a  preliminary  comparison  with 
existing  experimental  data,  a  penetration  model  of  the  form  developed 
herein  shows  some  promise  for  predicting  impact  behavior  over  a  wide 
velocity  range  and  for  different  projectile  and  target  materials. 
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-  constants  (Eq.  21) 

-  constants 

-  series  of  constants  as  C^,  C2»  •  .  . 

-  dilatational  wave  velocity 

-  projectile  diameter 

-  energy 

-  impact  kinetic  energy  of  projectile 

-  tensile  elongation  (strain)  at  fracture  {%) 

-  target  resistance  force 

-  modulus  of  rigidity 

-  shear  strain  at  fracture  [%) 

-  Brinell  hardness 

-  bulk  modulus 

-  constant 

-  mass 

-  crater  depth  measured  from  original  target  surface 

-  pressure 

-  radius 

-  shear  strength 

-  target  temperature 

-  time 

-  ultimate  strength  (tensile) 

-  shock  velocity 
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-  crater  volume  below  original  target  surface 

-  velocity,  impact  velocity 

-  yield  point  or  yield  strength 

-  shear  strain 

-  Poisson's  ratio 

-  mass  density 

-  shear  stress 


B 

c 

i 

P 

pr. 

S 

t 

0 

1 

1,2, 3 •  •  • 


-  backsplash 

-  crater 

-  interface  between  projectile  and  target 

-  projectile 

-  pressure 

-  denotes  material  strain,  as  Eg  (strain  energy) 

-  target 

-  denotes  undisturbed  (ahead  of  shock) 

-  denotes  compressed  condition  (behind  shock) 

-  used  with  constants 
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INTRODUCTION 


The  purpose  of  this  study  is  to  gather  and  assemble  existing  data 
on  ballistic  impact  and  on  material  failure,  especially  at  high  impact 
velocities  or  large  loading  -  to  establish  the  relative  importance  of 
such  factors  as  projectile  velocity,  mass,  projectile-target  strengths, 
ductilities,  densities,  compressibilities,  etc.,  and  to  use  this  infor¬ 
mation  to  deduce  the  mathematical  relationships  of  critical  factors  as 
the  target  structure  responds  to  impact  and  is  penetrated. 

Existing  experimental  data  relative  to  ballistic  impact  at  high 
velocities  are  being  evaluated  on  a  statistical  bases  through  the  use 
of  an  RPC  4000  digital  computer.  The  general  form  of  the  statistical 
approach  was  outlined  in  the  First  Quarterly  Report  ^ .  A  preliminary 
analysis  of  the  correlation  between  depth  of  penetration  and  ten  inde¬ 
pendent  variables  was  reported  in  the  Second  Quarterly  Report  2  and 
discussed  more  fully  in  the  Third  CXjarterly  Report  3. 

In  addition,  the  general  areas  of  ballistic  impact  and  material 
failure  are  being  investigated  in  order  to  develop  relationships  which 
may  be  tested  against  existing  experimental  data.  Some  general  aspects 
of  target  behavior  under  ballistic  impact  were  discussed  in  the  Second 
Quarterly  Report,  and  some  justification  was  given  for  the  use  of  static 
or  quasi-static  material  parameters  in  the  initial  statistical  analysis. 
Initial  attempts  to  formulate  a  theoretical  model  for  the  purpose  of 
testing  accumulated  experimental  data  were  outlined  in  the  Third  Quarterly 
Report.  The  basis  of  the  theoretical  work  was  the  energy  conversion  pro¬ 
cesses  occurring  during  impact. 


EMPIRICAL  MODEL 


An  analysis  of  the  penetration  depth  in  semi-infinite  targets  divided 
by  the  diameter  of  the  projectile  (Pc/Dp)  has  been  completed  on  1,270  shots 
which  were  reported  in  enough  detail  to  assign  strength  parameters  to  the 
target  materials.  The  experimental  data  was  divided  into  two  parts: 

(l)  those  having  impact  velocities  below  V3  (Hopkins  and  Kolsky  4),  and  (2) 
those  with  impact  velocities  near  or  above  Vg,  where 


The  distribution  of  target  materials  as  a  function  of  impact  velocity  is 
given  in  Figure  1  for  the  low  velocity  set  (985  shots)  and  in  Tigure  2  for 
the  high  velocity  set  (287  shots). 
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Figure  1.  Distribution  of  Target  Materials  as  a  Function 
of  Velocity  -  985  Low  Velocity  Shot* 
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Figure  2.  Distribution  of  Target  Materials  as  a  Function  of  Velocity 
297  High  Velocity  Shots 
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As  was  done  in  the  previous  analysis  of  (Pc)>  reported  in  the  Third 
Quarterly  Report  3,  (pc/^p)  was  fitted  to  a  simple  power  law  formula 


with  no  additional  assumptions  being  made.  This  equation  was  then  reduced 
to  its  associated  linear  form 


^  Vc/Dr  ~  Lo^-k0  i~  h  f  Ley  IT  +-  C<rj.  /Ot  +- . 

and  the  method  of  least  squares  was  used  to  determine  a  "best”  set  of  values 
for  the  variable  exponents  (coefficients)  k^  based  on  minimizing  the  value 
of 

£  ^  )c alc Ja4.J  3 

The  independent  variables  were  dropped  one-by-one,  starting  with  the  last, 
and  the  effect  on  the  remaining  kj  determined.  These  are  summarized  in 
Tables  I,  II,  III  and  IV  for  the  low  velocity  data  and  Tables  VI,  VII,  VIII, 

IX,  X,  and  XI  for  the  high  velocity  data. 

After  the  k.  were  computed,  the  usual  statistical  tests  were  made  to 
determine  the  validity  of  the  results.  The  per  cent  of  variance  in  log 
Pc/Dp  that  is  explained  by  the  variation  in  the  logs  of  the  independent 
variables  is  given  by  the  multiple  correlation  coefficients  listed  in  each 
table.  Individual  correlation  coefficients  between  all  variables  were  also 
computed  on  the  low  velocity  shots.  These  results  appear  as  Table  V. 

From  Table  V,  it  may  be  seen  that  the  close  correlation  between  the 
dilatational  wave  velocity  (c-t),  the  ultimate  strength  (U^),  the  shear  strength 
(S^),  the  Brinell  hardness  (H^),  and  the  yield  strength  (Y-j-)  render  these  in¬ 
dependent  variables  statistically  indistinguishable.  In  addition,  the  correla¬ 
tion  coefficient  between  the  bulk  modulus  (K^)  and  each  of  the  above  variables 
is  greater  than  0.50.  Any  of  the  empirical  power  law  equations  containing 
more  than  one  of  these  closely  correlated  variables  is  likely  to  be  completely 
misleading  in  the  coefficients  (exponents)  of  the  correlated  variables.  All 
the  rest  of  the  variables  including  the  per  cent  elongation  (e^-)  were  sta¬ 
tistically  independent. 

Although  the  results  of  the  above  analysis  have  not  been  thoroughly 
scrutinized,  a  few  general  observations  may  be  made.  First,  note  the  start¬ 
ling  difference  in  the  multiple  correlation  of  the  low  velocity  and  the  high 
velocity  data.  The  best  correlation  of  the  high  velocity  data  (49.5%)  is 
little  more  than  half  the  best  correlation  of  the  low  velocity  data  (85%), 
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in  spite  of  the  fact  that  the  high  and  the  low  velocity  distributions 
are  not  as  widely  separated  as  one  would  like.  It  should  be  noted  that 
the  high  velocity  correlation  has  not  been  completely  checked  for  errors 
which  could  considerably  alter  the  correlation  coefficients. 

It  is  interesting  to  note  that  as  long  as  five  or  six  of  the  inde¬ 
pendent  variables  are  used  in  the  empirical  formula  with  the  low  velocity 
data,  (pc/°p)  or  (Pc)  3  appears  to  be  a  function  of  the  first  power  of 
the  impact  velocity.  This  linear  dependence  of  (Pc)  on  (v)  agrees  with 
the  discussion  of  the  simple  power  law 

Pc=£uri-  (1) 


by  Herrmann  and  Jones  5  in  which  they  find  that  the  experimental  data,  at 
least  for  aluminum,  fits  n  =  1  over  the  approximate  velocity  of  7500  - 
12,000  feet/second.  They  also  point  out  that  in  the  velocity  range  of 
10,000  -  20,000  feet/second  penetration  may  be  approximated  by  Eq.  (l)  with 
n  -  2/3,  and  from  30,000  -  200,000  it  may  be  approximated  by  Eq.  (l)  with 
n  =  1/3.  Tables  VI  -  X  indicate  that  the  high  velocity  penetration  data 
used  in  the  present  empirical  analysis  is  indeed  best  fit  by  a  power  of 
velocity  slightly  less  than  2/3.  For  some  unexplained  reason,  Table  XI 
seems  to  indicate  a  dependence  on  a  power  of  velocity  nearer  to  1/3  than 
2/3. 


The  low  velocity  data  also  seem  to  indicate  a  dependence  of  (Pc/Dp) 
on  the  first  power  of  the  projectile  density  (/O  ).  This  coupled  witn 
the  linear  dependence  on  velocity  mentioned  abovepfits  Herrmann  and  Jones 
assumption  of  a  resistive  force  dependent  on  velocity  as 

F  *  v* 

which  leads  to  a  penetration  expression  of  the  form 


(2) 


where  k  is  a  constant.  As  they  point  out,  an  equation  somewhat  like  Eq.(2) 
can  be  used  to  fit  low-velocity  impact  data. 

The  high-velocity  data  show  a  dependence  of  ( Pc/Dp )  on  ( /op)  to  a 
power  slightly  greater  than  l/3.  This  coupled  with  the  dependence  on  the 
2/3  power  of  velocity  is  in  agreement  with  the  widely  used  empirical 
expression 


16 


(3) 


I 


J/3  2/3 

Pf  V 


which  is  obtained  for  the  case  where  the  crater  volume  is  proportioned  to 
the  energy  and  the  craters  are  hemispherical  in  shape.  The  fact  that 
Palmer  ^  finds  the  Utah  data  fits  an  equation  similar  to  Eq.  (3),  but  with 
the  exponent  of  ( /Op)  as  £,  is  also  understandable  on  the  basis  of  the 
high  velocity  data,  especially  in  Table  XI. 


Both  the  high  and  low  velocity  data  indicate  a  real  dependence  on  at 
least  one  target  strength  parameter  (c^.  and  also  being  classed  as  strength 
parameters  because  of  their  close  correlation  with  the  strength  parameters). 
The  low  velocity  data  seems  to  favor  a  dependence  of  penetration  on  Brinell 
hardness  to  the  -1/3  power  (Table  II),  which  agrees  with  Herrmann  and  Jones 
and  Palmer  6,  or  on  shear  strength  or  bulk  modulus  to  a  power  of  about  -J- 
(Tables  III  and  IV).  The  high  velocity  data  (Tables  VII  -  XI)  shows  a  depen¬ 
dence  of  penetration  on  Brinell  hardness,  shear  strength,  yield  strength,  or 
ultimate  strength  to  a  power  of  about  -0.2,  or  on  the  bulk  modulus  to  a  power 
of  about 


The  role  of  the  target  density  is  clearly  not  well  defined.  A  possible 
reason  for  this  is  given  in  the  discussion  of  the  theoretical  model.  The 
target  temperature  coefficients  are  not  considered  particularly  significant 
since  the  shots  in  which  temperature  was  varied  represent  such  a  small  per¬ 
centage  of  the  total. 


THEORETICAL  MODEL 


GENERAL 

The  study  of  the  transformation  processes  involved  in  the  dissipation 
of  Impact  kinetic  energy  in  a  hypervelocity  ballistic  impact  is  being  con¬ 
tinued.  A  semi-rational  penetration  expression  has  been  developed  from  a 
work-energy  consideration  which  suggests  that  the  nonrecoverable  target 
compression  and  shear  strain  energies  may  account  for  most  of  the  kinetic 
energy  of  the  projectile.  This  approach  is  not  intended  to  delineate  the 
relative  importance  of  the  mechanisms  of  cavitation,  shear,  plastic  flow, 
etc.,  but  rather  to  Indicate  that  the  resulting  cavity  formation  is  a 
function  of  certain  resistance  parameters  which  appear  to  be  primarily 
density,  compressibility  and  shear  toughness  of  th^ target  material  for  a 
given  set  of  impact  conditions.  This  estimate  of  the  situation  is  intended 
to  refer  primarily  to  the  more  ductile  targets  (as  most  of  the  experimental 
targets  have  been  -  -  aluminum,  steel,  copper,  lead,  etc.),  although  material 
shear  toughness  does  provide  a  measure  of  relative  ductility  or  brittleness. 
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The  kinetic  energy  of  a  projectile  can  be  considered  as  its  capacity 
to  do  work  as  a  result  of  its  velocity.  From  the  principle  of  work  and 
kinetic  energy,  it  can  be  said  that  the  total  work  done  by  a  projectile 
on  a  target  is  equal  to  the  change  in  kinetic  energy  of  the  projectile, 
since  the  negative  work  done  by  the  target  on  the  projectile  should  be 
equal  numerically  to  the  positive  work  done  by  the  projectile  on  the 
target.  Thus,  the  impact  kinetic  energy  should  be  equal  to  the  net 
positive  work  done  on  the  target. 

The  actual  work  done  in  moving  the  crater  material  backwards  is 
negligible  compared  with  the  work  required  to  fracture  or  dissociate  this 
material  from  the  target.  Any  recrystallization  that  takes  place  is  a 
result  of  internal  deformation  accompanied  by  heat.  It  appears  that  the 
net  work  done  on  the  target  may  be  determined  largely  by  the  nonrecoverable 
or  permanent  deformational  energy  absorbed  by  the  target  material.  That  is 
to  say;  the  energies  dissipated  as  heat,  sound,  electromagnetic  radiation, 
etc.,  are  essentially  a  result  of  material  deformation  (accompanied  by  frie-. 
tlonal  effects,  etc.)  and  that  the  permanent  deformation  itself  constitutes 
the  primary  work  done  by  the  projectile.  It  isreallzed  that  many  uncertain¬ 
ties  such  as  energy  going  into  shock-wave  formation  and  crack  propagation 
in  brittle  materials,  extent  of  inelastic  deformation  into  the  target,  etc,,* 
are  Involved,  although  these  features  are,  at  least  in  part,  a  function  of 
the  target's  ability  to  absorb  strain  energy.  Of  course,  there  are  bound 
to  be  other  considerations,  such  as  movement  of  the  entire  target,  etc., 
which  are  assumed  to  be  secondary  in  this  treatment. 

The  material  deformation  itself  seems  to  refer  to  the  primary  work 
being  done  on  the  target,  and  the  corresponding  deformational  energy 
absorbed  is  dissipated  through  the  different  processes.  It  is  thus  be¬ 
lieved  that  the  correlation  of  impact  kinetic  energy  with  either  mechanical 
deformational  work  or  the  combined  mechanical,  chemical,  acoustical,  thermal, 
electrical,  etc.,  energies  of  dissipation  provides  dual  methods  for  analyzing 
the  impact  problem.  In  the  latter  approach  the  target  strain  energy  absorbed 
and  any  work  or  energy  required  to  recrystallize  part  of  the  target  should 
not  be  combined  with  the  sum  total  of  the  backsplash,  heat,  sound,  radiation, 
etc.,  energies  eventually  lost  from  the  target.  In  connection  with  the 
recrystallization  process,  the  work  or  energy  required  to  recrystallize  the 
grains  might  be  considered  a  part  of  the  deformational  work  being  done  on 
the  target.  Presumably,  this  effect  could  be  included  in  the  general  defor¬ 
mational  parameters  chosen. 

The  problem  of  bridging  the  gap  between  lower  velocity  and  hypervelocity 
ballistic  impact  theory  requires  an  explanation  of  how  the  projectile  moving 
into  a  target  with  an  initial  velocity  greater  than  that  of  any  stress  wive 
that  can  be  detached  affects  the  cratering  process.  One  point  that  seems 
evident,  however,  i6  that  in  both  cases,  the  cratering  process  depends  largely 
upon  the  target's  ability  to  absorb  energy  which  must  be  related  in  some  way 
to  pertinent  impact  and  material  parameters.  If  the  relations  governing  the 
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pressures  developed  In  subsonic  and  hypervelocity  impacts  constitute  the 
primary  difference  in  the  two  cases,  it  is  shown  in  the  discussion  of 
Eqs.  (10,  (19),  (27),  and  (63)  that  the  subsonic  pressures,  determined 
from  the  Bernoulli  relation  (Eq.(lO))  and  the  hypervelocity  shock  pressures, 
determined  from  one-dimensional  shock  theory  (Eq.  (19)),  yield  similar  pene¬ 
tration  expressions  in  all  but  possibly  extremely  high  velocity  impacts. 

It  is  suggested  that,  in  general,  parameters  which  provide  a  measure 
of  a  target's  ability  to  "stop"  a  projectile  (and  in  doing  so  to  describe 
the  cratering  process)  are  appropriate  compression  and  shear  energy  ab¬ 
sorbing  parameters.  The  energy  absorbed  in  direct  compression  and  the 
shear  energy  absorbed  primarily  in  the  cratering  process  no  doubt  take 
on  varying  degrees  of  importance,  depending  on  target  material,  impact 
velocity,  etc.  There  are,  of  course,  discrepancies  involved  in  the  use 
of  static  or  low  strain-rate  material  properties,  but  these  difficulties 
are  presently  unavoidable  regardless  of  the  method  of  analysis.  It  does 
appear  that  the  pertinent  static  material  properties  are,  generally  speak¬ 
ing,  indicative  of  the  behavior  under  conditions  of  high  pressure,  high 
strain  rates,  etc.  That  if,  higher  pressures  generally  produce  greater 
strengths  and  ductilities  (at  least  shear  ductility)  but  less  compressibil¬ 
ity,  etc.,  and  higher  strain  rates  generally  produce  greater  strengths, 
but  less  compressibility  and  ductility,  etc. 

In  order  to  study  the  behavior  of  the  target  during  the  impact  process, 
it  is  necessary  to  consider  the  nature  of  the  loading  and  the  state  of  the 
material.  The  ordinary  use  of  the  terms  solid  and  fluid  seem  to  be  insuf¬ 
ficient  in  describing  the  target  material  under  such  high  strain  rates  and 
pressures,  etc.  The  hydrodynamic  theory  assumes  that  the  target  material 
has  been  stressed  far  beyond  its  strength  and  therefore  behaves  as  a  liquid 
with  virtually  no  shear  strength.  However,  since  extremely  high  pressures 
tend  to  congeal  liquids,  the  high  pressures  produced  during  impact  may  be 
siaqply  a  result  of  the  quasi-solid  target  material's  ability  to  develop 
such  high  pressures.  Also,  high  pressures  are  known  to  substantially  in¬ 
crease  shear  strengths  so  that,  even  though  some  type  of  material  flow  un¬ 
doubtedly  takes  place,  the  shear  strain  energy  absorbed  may  be  a  significant 
factor.  Regardless  of  the  magnitude  of  the  impact  velocity,  there  must  be 
a  certain  period  at  the  end  of  the  crater  formation  period  in  which  material 
strength  plays  a  dominant  role. 

The  compression  and  shear  strain  energy  absorbed  in  the  Impact  process 
depend  generally  on  the  following  factors: 

1.  The  magnitude  and  distribution  of  the  pressures  that 

the  target  material  develops  under  the  impact^conditions  along  with 
the  rate  at  which  the  high  initial  pressures  are  attenuated. 

2.  The  compressibility  of  the  target  material  or  its 
reciprocal,  the  ability  of  the  material  to  resist  compressive 
deformation. 
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* 


3.  The  shear  toughness  of  the  target  material  under 
the  impact  conditions. 

DEVELOPMENT  OF  MODEL 

The  general  expressions  for  compression  and  shear  strain  energy  ab¬ 
sorbed  for  the  case  of  a  simplified  hemispherical  surface  are* 

Compression  Strain  Energy 


(4) 

Shear.. Strain  Energy 


(5) 
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As  soon  as  the  impacting  materials  begin  to  deform,  the  maximum 
developable  pressure  might  be  considered  as  the  dynamic  fluid  pressure 
(or  stagnation  pressure)  for  the  smaller  density  ataterlal. 


f  e  c.  C/ov-2-; 


(6) 


where  C,  is  a  constant  and  is  the  smaller  of  the  projectile  or  target 
material  density. 

However,  a  better  estimate  of  the  pressure  developed  might  be  obtained 
for  the  case  in  which  both  projectile  and  target  are  regarded  as  incom¬ 
pressible  fluids.  In  this  case,  from  the  Bernoulli  relation, 

1°  =  z  /°t  V*  =  £  /0f  (V  -  irc  f  (7) 

where  v  =  impact  velocity 

vi  =  velocity  of  projectile-target  Interface 

From  Eq.  (7), 


(8) 

or 

V-  = 

(9) 

Hence, 

or 

f  'r/5^1  -  C  — —  /°fc -v-  z 

(10) 

II 

aP 

> 

1 

N 

(11) 

where  C2  is  a  constant  and  refers  to  a  composite  projectile- 

density  term  defined  as  ' 

■target 

^  =  Y°tY* 

(12) 
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Interesting  comparisons  are  made  on  Lines  5  and  6  of  Table  XII  which  show 
the  variation  of  the  composite  density  term  with  the  projectile  and  target 
densities.  Note  that  when  the  projectile  and  target  densities  are  equal, 
the  composite  density  term  is  0.25  times  either  one. 

Alternately,  the  initial  contact  pressure  might  be  determined  by  one- 
dimensional  shock  theory.  The  conservation  of  mass  and  momentum  equations 
are  written  for  a  stationary  shock  front,  where  the  undisturbed  condition 
(ahead  of  the  shock)  and  the  compressed  condition  (behind  the  shock)  are 
denoted  by  the  subscripts  0  and  1,  respectively. 


Conservation  of  Mass  yO,  ( U  -  V,  )  -  ^/0o  (  LL  —  Y0  ) 

(13) 

Conservation  of  Momentum  —  jpo  )-/o0  (u.-v0)(y,  -  Vo) 


where  u  is  the  shock  velocity  and  v  is  the  particle  velocity. 
Solving  Eq.  (13)  yields 


u.  =  v.  ±  M 

/ (/*■  -p* 


(u) 


and 


-y,  =  Vo  t 


(f>.  ~  Po)  (/°,  ~/0o  ) 


(15) 
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Table  XII.  olictlliMwt  CalcuUtltni 
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Hater lal  - 

4140  Stool 

Copper 

7|nc 

•aqn. 

Aiwa. 

Lead 

Reference 

1.  f 

raa  “boat  fit"  curve*  (Utah,  p.14)  — 

Er 

-10 

0.70  .  10 
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If  the  chang*  in  compressibility  with  pressure  is  neglected  (equivalent  to 
neglecting  the  nonlinear  ter*  in  Eq.  (21)),  the  shock  velocity  u  can  be 
taken  as  the  dilatational  wave  velocity,  or 


(16) 


The  appropriate  expression  for  the  dilatational  wave  velocity  c  for  the 
case  where  the  compressibility  behavior  is  assumed  to  be  linear  and  where 
the  contribution  of  shear  stiffness  G  is  neglected  (inelastic  range)  is 


Solving  Eqs.  (14),  (15)  and  (16),  following  the  assumption  that  VeK0 
and  -p0=o  (corresponding  to  the  initial  impact  condition)  yields, 


for  the  target 


~  “Vf  C±  /°-t 


(17) 


and  for  the  projectile 


f,  =  (  V  -  V,  )  C?/Or 


(18) 


where  unsubscripted  v  is  impact  velocity.  Eliminating  vj  between  Eqs. 
(17)  and  (18)  yields  an  expression  for  the  initial  interface  pressure 


f>,  = 


cv/°r  +  c*/0* 


(19) 


A  comparison  of  the  Bernoulli  pressure  (Eq.  (10))  with  the  shock  pressure 
(Eq.  (19))  in  the  Discussion  Section  shows  that  similar  penetration  ex¬ 
pressions  are  obtained  in  all  but  possibly  extremely  high  velocity  impacts. 

The  factors  that  govern  the  rate  at  which  the  initial  pressures  are 
attenuated  are  not  known.  Also,  the  distribution  of  pressure  around  the 
crater  rim  at  a  given  instant  is  unknown.  However,  a  representative 
pressure  variable  might  be  considered  to  be  a  function  of  the  initial 
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pressure  p.  Although  the  extent  of  the  Inelastic  deformation  in  the 
target  is  uncertain,  it  would  appear  to  be  generally  a  function  of  the 
target’s  ability  to  absorb  compression  energy.  Thus,  the  upper  limit 
for  the  integral  (Eq.  (4))  might  be  considered  as  some  function  of  the 
crater  depth  (constant  times  crater  depth  P_  for  linear  function). 
Assuming  this  function  to  be  linear,  a  solution  for  Eq.  (4)  would  be 


fesL  =  tv 


where  C3  is  a  constant. 

aV 

A  measure  of  the  unit  volume  change  of  the  target  material,  y  , 
might  be  taken  from  hydrostatic  compressibility  data  which  has  been  ob¬ 
tained  for  a  few  materials.  An  expression  of  the  form 


Af»-  V 


seems  to  fit  these  curves  fairly  well,  where  A  and  B  are  constants.  Actually, 
the  hugoniot  relationship  developed  for  a  given  equation  of  state  could  be 
used  here,  although  this  would  hardly  seem  consistent  with  the  use  of  static 
or  low  strain-rate  strength  parameters  (especially  when  the  strength  effect 
is  considerable). 

It  can  be  noted  in  Table  XII  that  the  factor  A  in  Eq.  (21 )  and  the 
compressibility  l/K  (reciprocal  of  bulk  modulus)  are  almost  identical  and 
that  the  factor  B  varies  similarly  between  materials.  If  the  non-linear 
term  in  Eq.  (21 )  is  discarded,  this  equation  can  be  taken  as 

AV  _  _L  ^ 

—  -  Kf  (22) 

Considering  the  shear  toughness  of  the  crater  material  to  be  the 
product  of  S  and  g  and  the  plastic  shear  strain  in  the  vicinity  of  the 
final  crater  boundary  to  be  a  function  of  this  toughness  factor;  also  assum¬ 
ing  the  extent  of  the  inelastic  shear  strain  into  the  target  to  be  some 
linear  function  of  the  crater  depth,  a  solution  of  Eq.  (5)  is 


c, 


where  C4  is  a  constant. 


25 


Equating  the  impact  kinatic  eaeray  to  tha  sub  of  the  two  work  or 
anargy  expressions,  Eqs.  (20)  and  (23),  yields 

Er  =  i  rrVr  -  T?’[c5  Afl  +  C,(S3)J  (24) 


(a) 

where  C5  and  C6  are  new  constants,  since  Vc  =  ^jTr‘^3  for  the  hemispheri¬ 
cal  crater.  Substituting  Eqs.  ( 11)  and  (22)  into  Eq.  (25)  yields 


which  will  be  used  to  investigate  the  Utah  data 
Eqs.  (11)  and  (22)  into  Eq.  (24)  yields 


(26) 


7 


Also  substituting 


The  use  of  Eq.  (21)  instead  of  Eq.  (22)  in  obtaining  Eqs.  (26)  and  (27) 
leads,  respectively,  to 


I 

Or  ~  L  +  Ck( 


(28) 
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and 


-  C7Bt/o»t'lr,'  +  C8(S3)tV^ 


(29) 


7 

The  Utah  data  will  now  be  used  to  study  the  proposed  penetration 
equation  in  the  form  of  Eq.  (26). 


=  - - -  (26) 

er  +  o-l  (sa)t 

This  form  of  the  equation  involving  the  simpler  compression  energy  ex¬ 
pression  (employing  Eq.  (22)  instead  of  Eq.  (21))  is  used  in  estimating 
the  relative  effects  of  compressibility  and  shear  toughness.  The  various 
values  are  shown  in  Table  XII  along  with  appropriate  reference  comments. 
One  shortcoming  (noted  on  Line  13  of  Table  XII)  is  the  necessity  of  using 
tensile  elongation  values  instead  of  ultimate  shear  strains  g  .  This  is 
not  too  unreasonable,  however,  since  both  provide  a  relative  measure  of 
ductility.  The  ultimate  shear  strains  for  the  different  target  materials 
could  not  be  found. 

The  Utah  report  7  showed  an  approximate  correlation  of  volume/energy 
with  shear  strength  St  (Eq.  (30))  for  all  target  materials  tested  except 
lead. 


X  ^  I 

FjO  Const.  (S)± 


(30) 


These  constants  are  repeated  on  Line  17  of  Table  XII  where  maximum  devia¬ 
tion  from  the  mean  value  is  70%.  Since  lead  is  the  most  ductile  material 
of  the  group,  the  possibility  of  a  closer  correlation  of  volume/energy 
with  shear  toughness  is  suggested  as 


X.  _  i 

E  f  Const.  ( 


(31) 


These  constants  are  shown  on  Line  18  of  Table  XII  where  the  maximum  de¬ 
viation  from  the  mean  value  is  only  22%,  using  the  tensile  elongation 
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parameters  for  the  ultimate  shear  strain  . 
Rewriting  Eq.  (26), 


a'  +  Ip  (ssV^  =  1  <32> 


or,,  referring  to  Table  XII, 

(  livne  |4  )  a,  -j-  (  l>ne.  IQ,)  O.j  =  I  (33) 


Approximate  values  of  a,  =  0.02  and  a,-  =  60  were  found  by  iteration  using 
the  Utah  data.  IVhcn  substituted  into'Tq.  (33),  these  constants  yield  a 
maximum  deviation  from  the  mean  of  20%  for  the  different  materials,  as 
shown  below. 


(34) 


aluminums  .406  a^  +  .0130  a0  —  0.8 

lead:  8.419  a-^  +  .0103  a3  ”=  1.1 

It  is  interesting  to  note  in  the  above  relations  that  the  very  much  higher 
coefficient  of  a-^  for  load  may  indicate  the  relative  importance  of  the 
density  term  when  the  impact  velocity  is  in  the  vicinity  of  the  target 
dilatational  wave  velocity  (the  case  only  for  lead  here).  The  maximum 
impact  velocities  were  about  2  km/sec. 

OTHER  MODELS 

This  section  is  adapted  from  the  excellent  paper  by  Herrmann  and 
Jones  5  for  purposes  of  comparison  and  discussion.  By  referring  to  an 
arbitrary  target  resistance  force  F,  the  work-kinetic  energy  expression 
can  be  written  as 

/Of>  1^,  r  dv  =  “F  d? c  (35) 


steel:  .074 
copper:  .322 
zinc:  .633 
magnesium:  .233 


aj  +  .0187  a-,  =  1.1 

a1  +  .0190  a.  =  l.l 
+  .0128  a-,  =  0.8 

a,  +  .  17b  a,  =  1.1 
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which  can  be  integrated  to  yield  a  penetration  expression  for  a  given 
definition  of  F.  This  approach,  however,  does  not  intuitively  include 
any  work  done  on  the  target  in  the  form  of  shear  deformation. 

Virtually  any  power  expression  of  the  form 


(36) 


can  be  obtained  from  Eq.  (35)  by  assuming  the  resistance  force  to  be  of 
the  form 


F  -  C, 


V 


J) 


(37) 


where  X  (  X  =  ~Z  yields  a  logarithmic  relation).  For  example,  if 

the  resistance  force  is  assumed  to  be  proportional  to  penetration,  as 

F  =  Crf  (38) 

or  proportional  to  all  four  parameters,  as 


F  =  C. 


Pc' 


Vr/°r 


or 


(39) 


Eq.  (35)  reduces  in  both  cases  to 


c  _ 


Y  V3 

f 


=  C 


nr 


(40) 


when  the  constant  of  integration  is  assumed  zero. 


If  C4  is  defined  as  some  target  material  strength  parameter,  one 
possibility  that  provides  an  approximate  data  fit  is 
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where  the  strength  parameter  R(.  is  the  Brinell  hardness.  Eq.  (40)  then 
becomes 


JL 

V  '<=> 


(42) 


(43) 


Eq.  (43)  can  also  be  obtained  by  equating  the  volume/energy  ratio  to 
and  assuming  a  hemispherical  crater. 


A  closer  fit  than  Eq.  (43)  has  been  found  by  adjusting  the  density 
ratio  power  as 


This  expression  involving  both  density  and  strength  terms  was  shown  by 
Herrman  and  Jones  to  provide  the  best  power  function  fit  for  a  variety  of 
target  and  projectile  materials  (a^  about  0.36),  although  still  over  a 
somewhat  limited  range  of  impact  velocities. 


If  now  the  resistance  force  F  is  assumed  to  be  dependent  on  the 
inertia  forces  exerted  by  the  target, 

d-  F  —  (  dm  CL  -  (^  Are*-  d.S 


=  Cc  Vf/2>  v  dt 

-  C  6  Vr^  /A>t  'V  dir 


dir 

dft 


(45) 
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and 


r=j :c»V4y»(/  (46) 

which  when  substituted  into  Eg,  (35)  yields  upon  integration 

V"^C7  S?J/n-V'  +  CB  <47> 

If  Cg  is  interpreted  as  a  material  strength  parameter,  independent  of 
velocity,  another  expression  similar  to  Eq.  (46)  might  be  written  as 


=  ClcHt) 


which  leads,  instead  of  Eq.  (47),  to 


Xs  z'c.  + 


Assuming  the  constant  cjj_  to  be  zero. 


v/3  aVA 


(%)  +  M 


Eq.  (50)  was  written  by  Herrmann  and  Jones  as 

ky  c"  ( a 


where  C^qH.  multiplier  (a  constant  for  a  given  material)  was  deleted. 
From  Eq.  (51),  a  least  squares  fit  of  existing  experimental  data  to  an 
equation  of  the  form 
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determined  a  more  appropriate  expression  than  Eq.  (51)  to  be 


which  can  be  approximated  by  the  expression 


for  a  velocity  range  limited  to  medium  and  moderately  high  impact  velo¬ 
cities.  Eq.  (54)  compares  very  closely  with  Eq.  (44). 

The -logarithmic  expression  involving  both  density  and  strength  terms 
has  been  shown  by  Herrmann  and  Jones  to  fit  experimental  data  over  the 
entire  experimental  range,  but  only  when  appropriate  values  for  the  con¬ 
stants  are  used  for  given  target  and  projectile  materials.  These  con¬ 
stants  were  found  to  vary  widely  between  materials. 

DISCUSSICN 

In  referring  to  theoretical,  semi-rational  or  entirely  empirical 
penetration  prediction  models,  one  point  that  must  remain  clear  is  that 
the  actual  mechanisms  involved  in  the  crater  formation  process  in  hyper- 
velocity  impacts  are  astronomically  complex.  Also  the  use  of  static  or 
relatively  low  strain-rate  (and  low  pressure)  mechanical  properties  of 
the  materials  renders  the  problem  that  much  more  vague.  It  can  only  be 
hoped  that  any  particular  penetration  expression  providing  an  engineering 
approximation  includes  those  factors  that  have  a  dominate  effect  on  the 
cratering  process  in  general. 

Specific  attention  will  now  be  focused  on  Eqs.  (26),  (27),  (44)  and 


where 


_  /°r  _ 

(  //y  +  J^t  ) 


In  addition  to  a  projectile  density  term,  Eqs.  (27),  (44)  and  (53)  all 
have  a  target  density  term,  a  target  strength  parameter  and  the  impact 
velocity  variable.  Whereas  the  Brinell  hardness  was  empirically  selected 
in  Eqs.  (44)  and  (53),  the  shear  toughness  was  rationally  included  in  Eq. 
(27).  The  form  of  the  density  terms  vary  widely  between  the  different 
expressions.  One  point  regarding  Eq.  (27)  is  that  the  sum  of  the  density 
(in  the  form  of  compressibility)  and  shear  toughness  terms,  each  with 
velocity  multipliers  provide  an  interesting  possibility  for  explaining 
behavior  over  wide  velocity  ranges.  For  example,  it  appears  that  the 
empirically  determined  constants  as  and  34  in  Eq.  (27)  would  require  the 
density  term  to  be  relatively  small  for  lvwcr  and  medium  velocities  when 
fitted  to  different  materials,  and  that  the  density  term  would  increase 
relative  to  the  toughness  term  as  the  velocity  increases  (note  the  two 
powers  of  the  velocity  variables  and  the  results  of^Eq.  (34)). 

A  linear  volume-energy  relationship  has  been  empirically  determined 
for  certain  ranges  of  impact  data. 
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(55) 


E  p  -  b  -h  A  VI 


or 

\v  __  ' 

FP  t/X  +  J  <56) 

where  b  refers  to  the  Ep  -intercept  and  d  is  the  slope. 

Eq.  (56)  has  the  same  form  as  Eq.  (26)  where  the  constant  slope  d 
in  Eq.  (56)  corresponds  to  the  constant  shear  toughness  term  in  Eq.  (26). 
Eq.  (56)  would,  however,  only  be  constant  when  the  b/Vc  term  was  negli¬ 
gible.  Actually,  the  Ep  -intercept  (or  b)  is  nearly  zero  for  most  lower 
and  medium  velocity  data.  This  then  adds  weight  to  the  form  of  Eq.  (26), 
or  its  equivalent  Eq.  (27),  where  the  density  term  is  small  for  lower  and 
medium  velocity  ranges.  For  higher  velocity,  the  constancy  of  Eq.  (56) 
would  no  longer  be  expected  and  would  no  longer  be  provided  by  Eq.  (26). 


From  an  analysis  of  Eq.  (27),  the  following  functions  are  written 
for  certain  limiting  conditions: 


JL  - 

Yr!'3 

T 


a  function 


of  (Sj) 


=  a  function  of 


=  a  function  of 


a  functi 


ion  of 


\ 

(57) 

jn 

I — 1 

i 

4- 

(58) 

+  EVP 

+  j_  _j_  +  J - ; 

+4yt  /V/fe  /Vfe  Pr0,‘ 

"1  (59) 

J  ) 
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a  function  of  ^  +/ {£'<■  - pf- -J  «*>> 


=  a  function  of 


(61) 


If  the  density  term  in  Eq.  (27)  is  negligible  at  lower  velocities, 
then  Eq.  (57)  should  be  approximately  correct.  The  penetration  relation 
with  the  shear  strength  to  the  minus  l/3rd  power  (Brinell  hardness 

gives  approximately  the  same  relation)  has  been  observed  in  the  past. 

It  is  shown  in  Table  XII  that  a  closer  penetration  relation  with  the  shear 
toughness  (Sg)-(-  results  for  a  certain  set  of  data  (actually  Vc/E 
proportional  to  (Sg)-j-  was  shown  in  Table  XII).  It  is  suggested  that  in 
empirical  correlations  the  product  of  and  gt  (the  shear  toughness) 
rather  than  the  two  parameters  separately  might  be  expected  to  provide 
bettor  data  fits. 


When  the  compressibility  (density)  term  is  negligible  in  Eq.  (27), 
Eq.  (58)  reduces  to 

Fwnc-f/on  of  V  ^  (62) 

In  low  and  medium  velocity  impacts  the  compressibility  effect  appears  to 
be  small,  and  the  2/3rds  power  has  been  found  to  govern  most  data.  It 
does  appear  that  in  these  velocity  ranges  the  shear  toughness  is  the  pre¬ 
dominate  parameter  compared  to  density  as  evidenced  for  example  by  a 
greater  penetration  in  lead  than  in  steel,  where  the  lead  shear  toughness 
is  smaller  than  steel,  although  the  density  is  larger. 


At  very  high  velocities  the  density  term  as  well  as  the  shear  term 
in  the  denominator  of  Eq.  (27)  would  presumably  be  significant,  so  that 
the  target  penetration  varies  functionally  with  the  velocity  as  Eq.  (58). 
This  would  tend  to  move  the  curve  of  Eq.  (62)  in  the  direction  of  the 
"hydrodynamic  theory"  l/3rd  power  relation. 


With  regard  to  the  question  of  the  appropriate  power  for  the  pro¬ 
jectile  density /Op  ,  values  of  /oJ/3  to  /O  %  or  so  have  been  reasonably 
well  established  as  limits  of  applicability.  For  Dqs.  (44)  and  (53) 


Herrmann  and  Jones  determined  values  of  approximately ,  while  the 
work  of  Palmer  ^  indicates  a  value  of pO  '/i  to  provide  the  best  fit. 

When  the  compressibility  term  is  negligible,  Eq.  (27)  yields  a  penetra¬ 
tion  expression  as  a  function  of ^Op'/-3  .  However,  as  the  compressibility 
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term  becomes  more  important!  Eq.  (59)  indicates  the  appropriate  /Op  -power 
to  move  in  the  direction  of  the  /O p ^3  -curve  (See  Figure  3). 

The  variation  of  penetration  with  target  density /Ot  is  shown  by 
Eq.  (60)  for  the  limiting  case  where  the  shear  toughness  term  is  small. 

The  first  term  is />t"^  which  is  then  adjusted  by  the  other  terms  in 
the  equation.  The  reason  a  /O^  “term  has  not  always  been  empirically 
included  in  the  past  could  be  that  the  density  term  is  relatively  small 
for  the  experimental  range  of  velocities!  and  thus  the  target  density 
effect  is  not  well  defined.  The  same  is  true  for  the  compressibility 
parameter  l/K^  given  by  Eq.  (61 ). 
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The  use  of  the  Bernoulli  pressure  given  by  Eq.  (10)  leads  to  Eq. 


(27) 
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[  asTT 


T  -f 


(S 


Vt 


J 


(27) 


If  now  the  shock  pressure  given  by  Eq.  (19)  is  used  instead  of  Eq.  (10), 
Eq.  (27)  becomes 


JL 

vv* 

r 


/3£ 


(63) 


One  apparent  difference  in  Eqs.  (27)  and  (63)  is  the  velocity  multiplier 
for  the  compressibility  term  only  in  the  case  of  Eq.  (27).  This  would 
simply  result  in  a  more  gradual  increase  in  the  compressibility  effect 
relative  to  the  shear  toughness  effect  as  impact  velocities  increase. 

Using  Eq.  (63)  instead  of  Eq.  (27),  Eqs.  (59)  and  (60)  become,  re¬ 
spectively 


64) 


-  a. 


F  o  nc.4-/  ovo  of 


(65) 


A  comparison  of  Eqs.  (64)  and  (65)  v/ith  Eqs.  (59)  and  (60)  shows  that 
very  similar  results  are  obtained  in  the  penetration  expressions,  as  far 
as  the  projectile  and  target  density  relationships  are  concerned.  Functions 
identical  to  Eq.  (65)  are  obtained  for  the  projectile  and  target  wave 
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velocities,  cp  and  c^.  It  follows  that  the  previous  statements  regarding 
the  relationship  between  penetration  and  target  density  apply  also  to 
projectile  and  target  sonic  velocities. 

Judging  from  the  previous  discussion,  an  equation  of  the  form  of 
Eqs.  (26)  or  (27)  seems  to  show  some  promise  as  a  penetration  expression 
for  predicting  impact  behavior  over  a  wide  velocity  range  and  for  different 
projectile  and  target  materials.  This  possibility  is  being  further  in¬ 
vestigated  statistically. 


FUTURE  WORK 


An  investigation  is  being  made  into  possible  causes  of  the  low  mul¬ 
tiple  correlation  coefficients  obtained  in  the  empirical  analysis  of  the 
high  velocity  data.  Discrepancies  in  the  data  can  be  spotted  by  calcu¬ 
lating  the  penetration  for  each  individual  shot  using  one  or  two  of  the 
better  empirical  equations  and  comparing  this  with  the  measured  penetra¬ 
tion.  In  this  way,  the  297  high  velocity  shots  used  in  the  present 
analysis  will  be  "cleaned  up".  Further  attempts  will  be  made  to  obtain 
new  high  velocity  data  and  thus  update  the  present  data. 

When  an  acceptable  set  of  high  velocity  shots  are  assembled,  this 
data  will  be  used  to  check  as  many  existing  theories  and  empirical  equa¬ 
tions  as  possible.  The  data  will  also  be  used  to  test  equations  developed 
under  the  present  contract. 

Further  theoretical  work  on  semi-infinite  targets  will  be  confined 
to  minor  adjustments  of  equations  already  derived  in  order  to  better  fit 
the  experimental  data.  It  appears  likely  that  the  present  program  will 
determine,  or  at  least  point  the  way,  to  an  acceptable  first  order  pene¬ 
tration  prediction  formula  for  semi-infinite  targets. 

Unfortunately,  the  contract  time  does  not  permit  an  analysis  of  thin 
targets  or  oblique  impacts.  It  is  hoped,  however,  that  this  work  may  be 
extended  in  order  to  update  the  present  analysis  as  new,  high  velocity 
data  appears  and  in  order  to  extend  the  study  to  thin  and  multiple  tar¬ 
gets  and  oblique  impacts. 
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